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NONLINEAR  EVOLUTION  OF  CONVECTING  PLASMA  ENHANCEMENTS  IN 
THE  AURORAL  IONOSPHERE  II:  SMALL  SCALE  IRREGULARITIES 

1.  INTRODUCTION 

Using  both  radar  and  satellite  measurements,  large  scale  convecting 
plasma  enhancements  in  the  auroral  ionosphere  have  recently  been  Identified 
and  studied  [Vickrey  et  al.,  1980].  Observed  in  regions  of  diffuse  auroral 
particle  precipitation  and  associated  field-aligned  currents,  these 

enhancements  have  overall  latitudinal  dimensions  of  a  few  hundred  kilometers, 
contain  relatively  steep  poleward  and  equatorward  edges,  and  have  been  shown 
to  be  approximately  field-aligned  resembling  vertical  slabs  of  ionization. 
Their  occurrence,  which  is  maximized  in  the  evening-midnight  sector,  is 
apparently  not  strongly  related  to  magnetic  activity  nor  to  E-region 

processes.  The  presence  of  plasma  density  irregularities  associated  with 
these  enhancements  has  been  verified  using  satellite  scintillation  studies 
[Fremouw  et  al .  1977;  Rino  et  al.,  1978;  Vickrey  et  al.,  1980].  The 
scintillation  data  have  indicated  that  the  electron  density  irregularities  are 
structured  like  L-shell  aligned  sheets  [Fremouw  et  al.,  1977;  Rino  et  al., 
1978].  In  addition,  Rino  and  Matthews  [1980]  have  shown  that  the 

scintillation  enhancements  resulting  from  these  irregularities  cannnot  be 
explained  in  terms  of  a  geometrical  enhancement  alone.  A  purely  geometrical 
enhancement  occurs  when  the  signal  propagation  path  intercepts  an  axis 
transverse  to  the  magnetic  field  along  which  axis  the  irregularities  have  a 
high  degree  of  spatial  coherence.  Moreover,  the  source  region  of  these 
scintillation  causing  irregularities  has  been  demonstrated  to  be  latitude 
limited  [Rino  and  Owen,  1980]  and  contained  in  a  vertical  slab  of  F  region 
plasma.  Using  simultaneous  rocket  probe,  scintillation  and  Incoherent 
scatter,  Kelley  et  al.  [1980]  have  also  recently  studied  several 

characteristics  of  auroral  F  region  irregularities,  e.g.,  total  electron 
content  and  spatial  power  spectra. 

Since  these  ionization  enhancements  have  been  observed  while  convecting 
equatorward,  their  poleward  edges  could  be  unstable  to  the  E.  x  gradient 
drift  instability  [Simon,  1963;  Li ns on  and  Workman,  1970]  as  observed  in 
artificial  ionospheric  plasma  clouds.  Indeed,  for  observed  [Vickrey  et  al., 
1980]  plasma  enhancement  density  gradient  scale  lengths  of  L  «  10-50  km  and 
convection  velocities  of  approximately  200  m/sec  (EQ  *  10  mV/m)  reasonable 
Manuscript  submitted  March  12,  1982. 
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growth  rates  for  the  E_  x  B_  gradient  drift  instability  can  be  expected 
since  y  *  a  (BL/cEQ)  *<  50-250  sec  where  y  is  the  E^  x  B_  growth  rate,  B  is  the 
ambient  magnetic  field  and  c  is  the  speed  of  light.  Moreover,  it  has  been 
shown  [Ossakow  and  Chaturvedi ,  1979]  that  by  applying  the  current  convective 
instability  [Lehnert ,  1958;  Kadomtsev  and  Nedospasov,  1960]  the  equatorward 
side  of  the  plasma  enhancements,  which  is  stable  to  the  _E  x  gradient  drift 
instability,  can  be  driven  unstable  by  the  ambient  field  aligned  particle 
precipitation  currents  in  conjunction  with  the  equatorward  density 
gradients.  Other  mechanisms  that  might  account  for  these  irregularities  are 
structured  low  energy  particle  precipitation  [Kelley  et  al.,  1980,  1982]  and 
irregular  field  aligned  currents.  Kesklnen  et  al.  [1980]  showed  that  the 
nonlinear  state  of  the  large  scale  irregularities  in  the  equatorward  edges  of 
these  plasma  enhancements  could  be  characterized  by  poleward  convecting  plasma 
depletions  and  equatorward-moving  enhancements.  In  addition,  they 

demonstrated  that  these  irregularities  could  be  described  by  inverse  power 
laws  in  the  nonlinear  regime.  Recently,  Keskinen  and  Ossakow  [1982]  discussed 
the  linear  stability  and  nonlinear  evolution  of  large  scale  convecting  plasma 
enhancements  in  arbitrary  ambient  electric  fields  in  the  auroral  ionosphere. 
These  studies  showed  that  convection  (through  the  E^  x  B^  gradient  drift 
instability)  is  the  primary  driver  of  long  wavelength  (3-100  km) 
irregularities  in  diffuse  auroral  F  region  plasma  enhancements.  However,  the 
aforementioned  satellite  scintillation  measurements  [Riuo  et  al.,  1978]  have 
Indicated  that  the  density  irregularities  associated  with  the  plasma 
enhancements  have  scale  sizes  down  to  hundreds  of  meters.  It  is  of  interest 
to  study  these  small  scale  (~  0.1-1  km)  irregularities  in  order  to  compare 
with  and  supplement  experimental  observations. 

In  this  report  we  present  a  linear  analytical  and  nonlinear  numerical 
study  of  small  scale  irregularities  applicable  to  local  unstable  regions  of 
large  scale  convecting  auroral  plasma  enhancements.  In  Section  2  we  give  a 
linear  stability  analysis  of  the  plasma  fluid  equations  which  describe  the 
evolution  of  density  fluctuations  in  the  auroral  F  region  ionospheric 
plasma.  In  Section  3  we  outline  the  methods  used  to  numerically  solve  these 
equations,  while  in  Section  4  our  principal  results  are  given.  Finally,  in 
Section  5  we  summarize  and  discuss  our  findings. 
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2.  EQUATIONS  OF  MOTION  AND  LINEAR  THEORY 


For  wavelengths  greater  than  the  ion  mean  free  path  we  use  fluid 
equations  to  describe  the  ion  and  electron  plasma.  The  following  geometry  is 
used:  the  y-axis  is  in  the  north-south  direction,  the  x-axls  points  west,  and 
the  z-axis  is  downward  along  the  magnetic  field.  In  this  report  we  ignore  the 
vertical  density  gradient  which  is  weaker  than  the  horizontal  plasma  density 
gradients  [Vickrey  et  al.,  1980]  in  the  typical  diffuse  auroral  plasma 
enhancements.  The  ion  and  electron  fluids  then  obey  the  following  equations 
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|f-  +  V*(n  v4)  *  0  (2) 
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Here  nfl  (o  ■  i  or  e)  is  the  species  density  and  E[  is  the  total  electric 
field.  Since  we  will  be  interested  in  low  frequency  fluctuations  we  have 
Ignored  inertial  terms  in  the  electron  and  ion  momentum  equations  (3)  and  (4). 
Equation  (5)  results  from  the  assumption  of  quasineutral 


fluctuations  n  *  n,  =  n.  In  addition,  vft  and  Vn  refer  to  the  electron  and 
e  i  u  u 

ion  velocities  along  the  magnetic  field  giving  rise  to  the  diffuse  auroral 

current*  The  symbol  v.  denotes  the  ion-neutral  collision  frequency,  v  .  the 

in 

electron-collision  frequency,  c  the  speed  of  light,  Tg  *  T  the  species 
temperature,  cg  the  ion  acoustic  speed  and  the  ion  (electron)  gyro- 

frequency.  We  have  neglected  vgn  compared  with  vfil  and  taken  «  1  for 

o  *  i,e  (F  region  approximation). 

Any  two  of  equations  (1),  (2),  and  (5)  provide  a  complete  description  of 
the  problem.  We  will  use  the  ion  continuity  equation  (1)  and  (5).  After 
separating  the  total  electric  field  into  an  ambient  and  fluctuating 
part  E  ■  -  V^6<J>  and  transforming  to  a  frame  drifting  with  velocity  -  - 

(c/B)  1 2.  x  Eg  -  < vln/S2i )  Eg]  we  can  write 
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where  *  t  (vQ  -  VQ).  Linearizing  (6)  and  (7)  by  separating  n  *  nQ(y) 
+  5n  with  6n,  6$  “  exp]i(kxx  +  kzz  -  lot)],  u>  ■  +  iy,  kL  »  1,  L  1  = 

(l/nQ)  (3nQ/3y)  we  find  a  growth  rate  (k(  =  k ;g) 
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where  0  =  kz/kx,  JEq  h  E^,  and  D±  -  <vei/fteCi)cs2  and  Dl  “  (cS2/vin^1  + 
f(vin/J5i)2/^(velvin/nefli)  +  (^z2/^i2))]}*  The  general  expression  for  the 
instability  growth  rate  y  using  arbitrary  directions  of  k^  and  Eq  can  be  found 
in  Reskinen  and  Ossakow  [1982] •  In  regions  of  plasma  enhancements  where 
3nQ/3y  <  0  (L  <  0)  we  find  the  condition  for  unstable  growth  to  be 


fl  1L| 
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where  we  have  taken,  for  example,  the  currents  to  be  downward,  i.e.,  <  0. 
The  effects  of  the  field-aligned  currents  will  be  able  to  reduce  (0  <  0)  or 
enhance  (0  >  0)  the  _E  x  gradient  drift  instability  growth  rate.  However, 
when  3nQ/3y  >  0  (L  >  0)  the  condition  for  unstable  growth  can  be  satisfied  for 
large  enough  current  velocities 
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The  expression  for  the  growth  rate  y  in  equation  (8)  can  be  maximized  as  a 

function  of  0  *  k./k  ,  a  measure  of  field-alignment,  using  3y/30|„  .  *  0 
I  X  vJ“U 
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Using  typical  diffuse  auroral  F  region  parameters  v^/n^  a  10  » 

v  . /SI  “  10  ^  E  *  10  mV/m,  J,  ■  n  eV.  »  1  pA/m2,  B  ■  0.5G, 
el  e  ^2  ox  i  o 

n  a  105cm  this  gives  |0  I  *  10  ,  i.e.,  approximate  field  alignment, 

o  ® 

Inserting  these  parameters  into  eq.  (8)  with  L  »  20  km,  *  0.2  mVsec 

and  D  a  108m2/sec  we  find  that  the  fastest  growing  linear  modes  have  growth 

*  .1  - 

times  y  *  10z  8ec. 
max 
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3.  NONLINEAR  THEORY 


In  order  to  study  the  nonlinear  evolution  of  these  small  scale 
Irregularities,  we  must  resort  to  a  numerical  solution  of  the  nonlinear  set  of 
equations  (6)  and  (7)  due  to  their  complex  nature.  Equations  (6)  and  (7)  can 
be  written  in  dimensionless  form  by  introducing  the  following  scaled 
quantities  n  =  no/NQ,  6$  =  6$/BL,  x  **  x/L,  y  *  y/L,  z  -  z/L,  t  **  ct/L  as 
follows  (where  we  have  dropped  the  tilde  for  clarity) 
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with  ct,  i  -  1,  ....  9  dimensionless  constants  given  by  cj  *  vin/0j» 
c 2  *  (vin/01)(T1/eBL)  +  (vei/fle)(cs2/nicL),  c3  -  c82/vincL» 

c4  "  neni^eV  cs  “  Eox/B’  Cfi  "  Eoy/B»  c?  *  (Oi/^XVj/c ), 
c $  -  T/eBL,  c9  -  (‘VV'eV  c8* 

In  the  following  numerical  simulations  we  take  advantage  of  the  fact  that 
the  fastest  growing,  most  dangerous  modes  from  linear  theory  are  almost  field-* 
aligned,  i.e.,  k(/k^  «  1  where  kj(k^)  is  the  component  of  k^  parallel 
(perpendicular)  to  the  magnetic  field.  These  waves  are  of  most  interest  to  us 
and,  as  a  result,  we  solve  equations  (10)  and  (11)  in  a  plane  containing  these 
modes  which  is  nearly  perpendicular  to  the  magnetic  field  while  fixing  the 
value  of  k(/k  «  1.  A  similar  approach  has  been  adopted  In  numerical  studies 
of  drift-wave  [Lee  and  Okuda,  1976]  and  trapped-particle  {Hatsuda  and  Okuda, 
1976]  instabilities  in  laboratory  plasmas.  The  system  of  equations  (10)  and 
(11)  was  first  transformed  to  the  x'y'z'  coordinate  system  (as  shown  in  Fig. 
1)  by  a  simple  rotation  about  the  y-axis  by  the  angle  0  »  k _ /k  «  1  using 
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where  0  is  the  angle  for  maximum  linear  growth  rate  defined  by  eq.  (12)  for  a 

definite  set  of  parameters  v.  /fl. ,  cE  /BV,,  v  /SI  .  Since  0  «  1  this 

r  in  i  ox  d  ei  e 

transformation  can  be  written  3/3x  =  3/3x',  3/3z  =  0  3/3x',  3/3y  =  3/3y'  with 

3/3z^  a  0.  As  a  consequence  the  three  dimensional  problem  is  reduced  to  two- 
dimensions.  By  solving  equations  (10)  and  (11)  in  the  x'y'z'  coordinate 
system  a  small  but  finite  k^  is  effectively  introduced  into  the  model. 

Equations  (10)  and  (11)  were  then  solved  numerically  on  a  mesh  consisting 
of  64  grid  points  in  the  north-south  direction  (y-direction)  and  64  grid 
points  in  the  east-west  direction  (x-di rect ion)  with  constant  grid  spacing  of 
15  m.  As  a  result,  the  simulation  plane,  which  is  taken  to  be  essentially 
horizontal  at  an  altitude  of  350  km  in  the  diffuse  auroral  F  region,  has  a 

north-south  and  east-west  extent  of  960  m,  respectively.  The  field  aligned 

currents  are  taken  to  be  constant  in  space  and  time  over  the  grid.  The  plasma 
density  n  in  equation  (10)  was  advanced  in  time  using  a  multi -dimensional 
flux-corrected  variable  timestep  leapfrog-trapezoid  scheme  [Zalesak,  1979] 
which  is  second  order  in  time  and  fourth  order  in  space.  At  each  timestep  the 
self-consistent  electrostatic  potential  6$  of  the  plasma  enhancement  in  eq. 
(11)  was  determined  using  a  Chebychev  iterative  method  [McDonald,  1980]  with  a 
convergence  criterion  of  10“^.  Since  we  are  considering  a  small  local 
unstable  region  (960  m  by  960  m)  of  a  large  scale  plasma  enhancement  which  is 
several  hundred  kilometers  in  extent,  periodic  boundary  conditions  were 
imposed  both  in  the  east-west  and  north-south  directions. 

4.  RESULTS 

In  the  following  we  consider  the  linear  and  nonlinear  evolution  of  small 
scale  irregularities  in  plasma  enhancements  in  the  diffuse  auroral  F  region 
ionosphere  in  an  approximately  horizontal  plane  at  350  km  altitude  almost 
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perpendicular  to  the  magnetic  field.  We  take  the  following  typical  parameters 

(Vickrey  et  al.,  1980;  Schunk  and  Walker,  1973;  Banks  and  Kockarts,  1973]  L 

=  20  km,  v.  /fl,  =  2  x  10~4,  v  . /SI  =  2  x  10“4,  =  10  mV/m  T  =  T.  = 

in  i  ei  e  ox  e  1 

1000°K  and  J.  =  1  pA/ra2  (which  gives  a  current  velocity  of  Vj  =  60  m/sec 

—3  a 

with  N  »  1  x  105cm  ).  In  addition,  we  assume  that  the  diffuse  auroral 
o 

particle  precipitation  current  J  is  downward  (V^  <  0)  and  spatially  and 

temporally  uniform  over  the  entire  plasma  enhancement.  In  order  to  find  the 

location  and  magnitude  of  the  maximum  linear  growth  rates  to  be  expected  with 

this  set  of  parameters  we  first  compute  0^  =  k||/kx  as  given  in  eq.  (9) 

with  V ,  =  -  l V ,  I  =  -  60  m/sec.  This  gives  two  values  for  6  which  are 

0  =  1.4  x  1U  and  0  =  -  6. 5  x  10  .  Using  eq.  (9)  and  considering 

wavelengths  X  =  2ir/k  =  500  m  the  first  value  0+  gives  a  maximum  linear 
X  X  -2-1 

growth  rate  yjaax  =  1.0  x  10  sec  on  the  poleward  side  (Sn^/Sy  <  0)  with 

linearly  damped  perturbations  Ymax  *  -3.5  x  10  ^sec  *  on  the  equatorward  side 

(  3n  /3y  >  0).  The  second  value  0  gives  only  a  marginally  unstable  growth 
°  -5  -1 

rate  of  y  =  8  x  10  sec  on  the  equatorward  side  with  damped  fluctuations 
max  -4-1 

Y  =  -3.7  x  10  sec  on  the  poleward  side.  These  results  agree  with  the 
m3  x 

experimental  observations  [Vickrey  et  al.,  1980]  that  the  largest  linear 
growth  rates  occur  on  the  poleward  side  of  the  equatorward  convecting  plasma 
enhancements.  In  this  case  the  effect  of  the  field-aligned  currents  is  to 
enhance  the  x  B_  gradient-drift  instability  growth  rate  on  the  poleward 

side.  The  current  velocities  are  too  weak  for  the  cases  studied 
observationally  to  give  appreciable  growth  on  the  equatorward  side  of  the 
plasma  enhancements.  We  will  then  consider  the  evolution  of  modes  satisfying 

0+  -  k„/k  =  1.4  x  10~5. 

I  x 

A  slab  approximation,  with  initial  density  profile  n  ( y ')  =  N  [1  - 

3  o  o 

y'/L  +  e(x',y')],  Nq  =  105cm  ,  is  used  to  model  a  small  local  unstable  region 

on  the  poleward  side  of  a  zero  order  equatorward  convecting  large  scale  plasma 
enhancement  in  the  diffuse  auroral  F  region  ionosphere.  In  this  assumed 

profile,  e(x',y')  denotes  the  initial  perturbation  and  L  =  20,  30  km  the 
initial  plasma  enhancement  density  gradient  scale  length.  We  consider  two 

models  distinguished  by  the  initial  seed  perturbations  e(x',y').  In  Model  1, 
purely  random  white  noise-like  initial  conditions  are  used  with  e(x',y') 
having  a  root  mean  square  value  of  0.01.  In  Model  2,  a  two-dimensional 

monochromatic  perturbation  is  employed.  For  Models  1  and  2  we  take  EQX  * 

10  raV/m,  EQy  *  0.  We  now  drop  the  prime  notation  for  clarity. 
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Figures  2(a)-2(c)  give  the  evolution  of  the  isodensity  contour  plot  of 
small  scale  density  fluctuations  $n(x,y)/nQ  using  Model  1  with  L  =  20  km. 
Figure  2(a)  shows  the  purely  random  nature  of  the  initial  conditions  still 
persists  at  t  =  100  sec.  Figure  2(b)  gives  the  evolution  of  6n/nQ  at  t  *  800 
sec  where  some  north-south  elongation  and  steepening  have  occurred.  Note  that 
local  density  enhancements  (6n/nQ  >  0,  solid  contours)  are  convecting 
northward  (poleward)  while  local  depletions  (6n/nQ  <  0,  dashed  contours)  are 
convected  southward  (equatorward ) .  This  relative  movement  of  enhancements  and 
depletions  is  directly  analogous  to  the  classical  Rayleigh-Taylor  instability 
in  a  heavy  fluid  supported  by  a  lighter  fluid.  Doppler  radar  backscatter 
signatures  [Hanulse  et  al.,  1981)  of  10.5  m  irregularities  in  the  evening- 
midnight  auroral  F  region  ionosphere  also  indicate  a  similar  convection 
pattern  in  a  plane  nearly  perpendicular  to  he  magnetic  field.  In  this 
experiments  southward  irregularity  convection  was  observed  approximately  1000 
km  due  north  from  northward  looking  HF  radar  while  farther  north  at  a  range  of 
1200  km  the  irregularities  appear  convecting  northward.  Finally,  Fig.  2(c) 
details  the  density  fluctuations  in  the  nonlinear  regime  at  t  =  1000  sec. 
Further  elongation  and  steepening  is  evident.  Similar  density  contour 

development  was  also  observed  for  the  other  plasma  enhancement  density 
gradient  scale  length  used,  L  *=  30  km. 

Figure  3a-b  give  sample  one-dimensional  spatial  power  spectra  in  the 
nonlinear  regime  at  t  *  1000  sec  both  in  the  east-west  (P(kx))  and  north-south 
(P(ky))  directions,  respectively  for  Model  1.  These  power  spectra  are  defined 
as  follows 

P(kx)  -  /  dkyP(kx,ky) 

and 

P(k  )  -  /  dk  P(k  ,k  ) 
y  x  x  y 

where  P(kx,ky)  =  (I^Ly)  * [6n(kx,k^ )/Nq ]2  is  the  spectral  density, 

6n  =  n-NQ  with  nQ  the  peak  plasma  enhancement  density,  and  l^Ly  is  the  area  of 

the  numerical  simulation  plane.  For  both  cases  these  power  spectra  are  well- 

fitted  with  inverse  power  laws  P(k  )  «  k  ""n  and  P(k  )  »  k  ~n  with  index 

x  x  y  y 

n  *  2-3.  The  spectral  Indices  are  in  agreement  with  those  obtained  from 
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Isodensity  contours  of  6n/nQ  with  L  ■  20  km  at  (a)  t  -  100  sec,  (b) 

t  ■  800  sec,  (c)  t  ■  1000  sec  Model  1.  Solid  contours  denote  local 

enhancements  6n/n  >  0;  dashed  contours  denote  local 

o 

depletions  6n/nQ  <  0.  The  observer  is  looking  upward  along  the 
magnetic  field  lines. 
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recent  scintillation  studies  [Erukhimov  et  al.,  1981]  of  small 
scale  (<  1  km)  plasma  turbulence  in  the  auroral  ionosphere. 

Figures  4(a)-4(c)  illustrate  the  evolution  of  the  density  fluctuations 

dn/n  for  Model  2  with  L  =  20  km.  The  simulations  were  initialized  with  a 
o 

general  monochromatic  two-dimensional  perturbation  of  the  form  [Rognlien  and 
Weinstock ,  1974;  Chaturvedl  and  Ossakow,  1979] 

6n(x,y)/n  =  Aj  jsin  k  y  cos  k  x  +  A2  Qsin  2  k  y 

°  «  y  x  >  y 

with  k  =  k  =  2n/960  m.  A;  ;  =  2  x  10  ,  A2  o  =  2  x  10  ^  where 

x  y  >  i 

A1  l(A2  o)  Is  linearly  unstable  (damped).  Figure  4(a)  gives  an  isodensity 
»  > 

contour  plot  of  6n(x,y)/nQ  at  t  =  0  sec.  The  initial  perturbation  describes  a 
sequence  of  local  enhancements  (Sn/n^  >  0)  and  depletions  (dn/n^  <  0) 
arranged  in  a  checkerboard  fashion.  Figure  4(b)  shows  the  evolution  of  fin/n^ 
at  t  =  900  sec  where  elongation  and  steepening  are  again  seen.  Figure  4(c) 
gives  Sn/n^  at  t  =  1200  sec  in  the  nonlinear  regime  where  further  steepening 
and  north-south  elongation  can  be  observed.  As  in  the  Model  1  the  late  time 
evolution  of  the  small  scale  plasma  density  irregularities  on  the  poleward 
side  of  large-scale  equatorward  convecting  plasma  enhancements  can  be 
characterized  by  poleward  moving  steepened  local  enhancements  and  equatorward 
convecting  depletions.  Spatial  power  spectra  in  the  nonlinear  regime  similar 
to  Model  1  were  also  observed  in  Model  2. 

However,  satellite  scintillation  studies  [Fremouw  et  al.,  1977;  Kino  et^ 
al. ,  1978;  Rlno  and  Matthews,  1980]  dealing  with  equatorward  convecting  large- 
scale  plasma  enhancements  have  indicated  that  density  irregularities  with 
scale  sizes  (0.1-1  km)  are  primarily  L-shell  (east-west)  aligned,  i.e.,  have  a 
higher  degree  of  spatial  coherence  along  L-shells  than  in  the  north-south 
direction.  On  the  contrary,  the  present  and  a  previous  study  [Keskinen  and 
Ossakow,  1982]  have  shown  that  density  irregularities  in  both  the  long  (1-100 
km)  and  short  (0.1-1  km)  wavelength  regime  in  equatorward  convecting  auroral 
plasma  enhancements  are  primarily  north-south  aligned  which  is  consistent  with 
the  nonlinear  development  of  the  Rayleigh-Taylor-like  E_  x  B_  gradient  drift 
instability.  As  shown  in  Model  2,  nonlinear  mode  coupling  effects  [Chaturvedi 
and  Ossakow,  1979]  do  not  appear  to  be  sufficient  to  account  for  the  L-shell 
alignment.  Keskinen  and  Ossakow  [1982]  have  studied  the  nonlinear  evolution 
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of  equatorward  convecting  plasma  enhancements  which  initially  contain  only  a 

north-south  density  gradient.  The  east-west  density  gradient  in  the  plasma 

enhancements  is  very  weak  [J.  Vickrey,  private  communication].  These  plasma 

enhancements  were  shown  to  be  unstable  and  to  break  up  into  primary  north- 

south  aligned  finger-like  structures  which  themselves  contain  sharp  east-west 

and  north-south  density  gradients.  If  these  long  wavelength  primary 

irregularities  have  a  component  of  convection  in  the  east-west  direction  (from 

a  north-south  electric  field),  then  secondary  smaller  scale  approximately  L- 

shell  aligned  structures  could  grow  on  the  east-west  density  gradients  of  the 

primary  irregularities.  A  north-south  electric  field  is  usually  present  in 

the  evening  diffuse  auroral  F  region  ionosphere  [Banks  and  Doupnik,  1975; 

Vickrey  et  al«,  1980].  (Some  evidence  that  this  process  can  occur  was 

presented  in  Keskinen  and  Ossakow  [1982].)  However,  the  east-west  convection 

must  be  of  sufficient  magnitude  compared  with  the  north-south  convection  to 

prevent  velocity  shear  stabilization  [Perkins  and  Doles,  1975]  of  the  smaller 

scale  secondary  irregularities.  For  an  initial  density  variation  along  the  y- 

direction  with  scale  length  L,  the  approximate  stability  criterion 

is  E  /E  >  2/kL  where  k  is  the  perturbation  wavenumber  and  E^XE,!  the 
oy  ox  r  °y  ox 

component  of  the  _E  x  _B  convective  electric  field  parallel  (perpendicular)  to 
the  initial  density  gradient. 

Figures  5(a)-5(c)  illustrate  the  evolution  of  the  density  fluctuations 
6n(x,y)/no  for  the  above  secondary,  two-step  model  in  which  smaller  scale  size 
irregularities  can  grow  on  larger  striation-like  structures.  We  assume  that 
the  large  scale  structures  are  in  the  nonlinear  regime  [Keskinen  and  Ossakow, 
1982]  and  evolve  on  a  slower  time  scale  than  the  smaller  scale  size 
irregularities.  The  initial  density  profile  which  describes  the  east-west 
density  gradient  of  a  small  local  region  of  a  large  primary  striation  is 
taken,  for  simplicity,  to  be  of  the  approximate  form  nQ(x)  =  N0 [  1  -  x/L 
+  e(x,y)]  with  L  =  5  km  and  e(x,y)  derived  from  white-noise  random  initial 
conditions.  For  simplicity,  we  take  the  total  electric  field  to  be  northward 
(E,  =  E^)  of  magnitude  |E^|  =  10  mV/m.  All  other  parameters  remain  the 
same.  Figure  5(a)  shows  the  evolution  of  the  isodensity  contour  plot  at  t  * 
100  sec.  Figure  5(b)  gives  6n/nQ  at  t  -  325  sec  where  some  steepening  and 
elongation  in  the  east-west  direction  (L-shell  alignment)  has  occurred. 
Finally,  Fig.  5(c),  which  shows  further  steepening  at  t  *  400  sec  in  the 
nonlinear  regime,  clearly  illustrates  L-shell  aligned  structures.  The  spatial 
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power  spectra  of  these  density  Irregularities  are  similar  to  those  found  in 
Model  1*  To  be  sure,  this  model  is  very  crude.  We  present  it  only  as  a  basic 
illustration  of  the  two-step  irregularity  generation  mechanism.  More 
realistic  models  will  be  studied  in  a  future  paper. 


5.  SUMMARY  AND  DISCUSSION 

We  have  performed  analytical  and  numerical  simulation  studies  of  small 

scale  (0.1-1  km)  irregularities  in  local  unstable  regions  of  large  scale 

convecting  plasma  enhancements  in  the  diffuse  auroral  F  region  ionosphere.  We 

have  shown  that  small  scale  size  density  fluctuations  can  be  destabilized 

primarily  on  the  poleward  sides  of  equatorward  convecting  plasma  slabs  by  a 

combination  of  the  effects  of  convection  and  field  aligned  currents.  In  a 

plane  nearly  perpendicular  to  the  magnetic  field  these  simulations  indicate 

that  this  destabilization  leads  to  steepened  striation-like  structures 

(elongated  in  the  north-south  direction  for  equatorward  convection)  which  can 

form  and  cascade  from  kilometer  to  tens  of  meter  scale  sizes  on  the  order  of 

an  hour.  The  c-e-dimensional  spatial  power  spectra  of  the  density 

irregularities  in  the  north-south  P(ky)  «  k  and  east-west  P(kx) 

«  k^  directions  can  be  described  by  power  laws  with  n  =  2-3  for 

wavelengths  2x/k  ,  2ir/k  a  80-960  m.  In  addition,  we  show,  using  a  very 
x  y 

simple  model,  that  the  experimentally  observed  L-shell  (east-west)  aligned 
nature  of  small  scale  (<  1  km)  irregularities  in  equatorward  convecting  large 
scale  plasma  enhancements  might  arise  from  a  secondary,  two-step  process.  In 
this  theory,  nonlinear  long  wavelength  (~  100  km)  primary  striations  can 
create  sharp  east-west  density  gradients  on  which  shorter  scale  size  (<  1  km) 
irregularities  can  grow  if  the  north-south  electric  field  is  of  sufficient 
magnitude  compared  to  the  east-west  field. 

In  this  study  we  have  examined  the  quasi  two-dimensional  linear  and 
nonlinear  evolution  of  small  scale  size  (0.1-lkm)  irregularities  in  local 
unstable  regions  of  larger  scale  (several  hundreds  of  kilometers)  convecting 
plasma  enhancements  in  the  diffuse  auroral  F  region  ionosphere.  This  has  been 
accomplished  by  solving  the  plasma  fluid  equations  in  a  horizontal  plane 
approximately  perpendicular  to  the  magnetic  field.  The  observed  plasma 
enhancements  are  three  dimensional  [Vickrey  et  al.,  1980] .  However,  the 
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horizontal  gradients  are  much  steeper  than  the  vertical  density  gradients 
allowing  one  to  approximately  model  the  plasma  enhancements  by  vertical 
slabs.  In  addition,  we  have  not  included  a  full  spectrum  of  finite 
k  modes  in  these  simulations.  However,  since  the  modes  with  maximum  linear 
growth  rate  have  kj/k^  «  1,  the  important  structuring  processes  will  occur  in 
the  plane  nearly  perpendicular  to  the  magnetic  field. 

Finally,  we  note  that  we  have  not  addressed  other  related  topics,  e.g. , 
the  source  mechanism  of  the  plasma  enhancements  and  their  coupling  to  the  L- 
region.  These  topics  will  be  discussed  in  future  studies. 
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OICY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  A  SPACE  CO.,  INC. 

3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT  DEPT  52-12 
OICY  ATTN  W.L.  IMHOF  DEPT  52-12 
OICY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 
OICY  ATTN  J.B.  CLADIS  DEPT  52-12 

LOCKHEED  MISSILE  &  SPACE  CO.,  INC. 
HUNTSVILLE  RESEARCH  4  ENGR.  CTR. 

4800  BRADFORD  DRIVE 
HUNTSVILLE,  AL  35807 
ATTN  DALE  H.  DIVIS 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
?.C.  BOX  5837 
ORLANDO,  FL  32805 
OICY  ATTN  R.  HEFFNER 

N.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  KA  02173 

OICY  ATTN  DAVID  K.  TOWLE 
OICY  ATTN  P.  WALDRON 
OICY  ATTN  L.  LOl'GKLIN 
OICY  ATTN  D.  CLARK 

MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
OICY  ATTN  X.  HARRIS 
OICY  ATTN  J.  MOULE 
OICY  ATTN  GEORGE  MROZ 
OICY  ATTN  W.  OLSON 
OICY  ATTN  R.W.  HALPRIN 


OiCY  ATT X  TECHNICS-  LIBRARY  SERVICES 


MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 


J.'TA  E 

ARBOR ' 

CA  93101 

~  ICY 

ATTN 

P.  rlSCr.ER 

•.  ICY 

.ATTN 

W.F.  C.-vEVIER 

C  ICY 

ATTN 

STEVEN  L.  GUTSCHE 

0  1  CY 

ATTN 

D.  S.AFFSNFIELD 

OICY 

AT  i  N 

R.  BGGUSCH 

OICY 

ATTN 

R.  HENDRICK 

OICY 

ATTN 

RALPH  KILB 

OICY 

ATTN 

DAVE  SOWLE 

OICY 

ATTN 

F.  FAJEN 

0  i  CY 

ATTN 

M.  SCHEIBE 

C  1  CY 

ATTN 

CONRAD  L.  LONGMIRE 

OICY 

ATTN 

WARREN  A.  SCHLUETER 

MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01 730 

01 CY  ATTN  JOHN  MDRCaNSTERN 
u ICY  ATTN  G.  HARDING 
OICY  ATTN  C.E.  CALLAHAN 

MITRE  COR? 

WESTGATE  RESEARCH  P.-JUK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
OICY  ATTN  W.  HALL 
C1CY  ATTN  W.  FOSTER 

PACIFIC-SIERRA.  RESEARCH  CORP 
1456  CLOVERFIELD  BLVD. 

SANTA  MONICA,  CA  90404 
OICY  ATTN  E.C.  FIELD,  JR. 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  T J  THIS  ADDRESS) 

OICY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PH J70METAICS ,  INC. 

442  MARRETT  ROAi 

N  G .  ON ,  0.,4/j 

G ICY  ATTN  IRVING  L.  KOFSKY 


PHYSICAL  DYNAMICS ,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  4  D  ASSOCIATES 

P.O.  BOX  9695 

MARINA  DEL  REY,  CA  9029) 

OICY  ATTN  FORREST  GILMORE 
OICY  ATTN  BRYAN  GABBARD 
OICY  ATTN  MILL I AM  B.  WRIGHT,  JR. 
OICY  ATTN  ROBERT  F.  LELEVIER 
OICY  ATTN  WILLIAM  J.  KARZAS 
OICY  ATTN  K.  ORY 
OICY  ATTN  C.  MACDONALD 
OICY  ATTN  R.  TURCO 

RANT)  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
■  ICY  ATTN  CULLEN  CHAIN 
OICY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  RC.AD 
SUDBURY  FLA  01776 

OICY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
SO  WEST  ENT  AVENUE 
NEW  YORK,  NY  10023 

OICY  ATTN  VINCE  TRAPANI 


SCIENCE  APPLICATIONS,  INC. 
P.O.  BOX  2351 


LA  JOLLA ,  CA 

92038 

OICY 

ATTN 

LEWIS  M. 

-INSON 

OICY 

A 1 1 N 

DAN  1 1.  l.  A . 

HAMLIN 

OICY 

ATTN 

E.  ERIE! LA 

C  ICY 

ATTN 

E.A.  STRAKER 

OICY 

A  A  i 

CURTIS  A. 

SMITH 

OiCY 

ATTN 

JACK  MC DOUG ALL 

SCIENCE 

APPL 

L CATIONS , 

INC 

1710  GOCDRIDGE  DR. 
'MCLEAN,  VA  22102 
ATTN :  J .  COCKAYNE 


PHYSICAL  DYNAMICS,  INC. 
P.O.  BOX  3027 
BELLEVUE,  KA  98009 
•OICY  ATTN  E.J.  FREMOUW 
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SSI  INTERNATIONAL 
333  RAVENSKOOD  AVENUE 


:nlo 

PARK,  l 

CA  94025 

01CY 

ATTN 

DONALD  N’EILSON 

c  :cy 

ATTN 

ALAS  BURNS 

C1CY 

C.  SMITH 

1 1  Y 

ATTN 

L.L.  COsb 

'?  1  C'i 

AT”' . 

DAVID  A.  JOHNSON 

C  ICY 

ATT  \ 

WALTER  G.  CHESKUT 

01CY 

ATTN 

CHARLES  L.  KINO 

01CY 

ATTN 

WALTER  JAYE 

01CY 

ATTN 

M.  BARON 

C1CY 

ATTN 

RAY  L.  LEADABRAND 

0  ICY 

ATTN 

G.  CARPENTER 

01CY 

ATTN 

G,  PRICE 

OiCY 

ATTN 

J.  PETERSON 

0  ICY 

ATTN 

R.  HAKE,  JR. 

01CY 

ATTN 

V.  GONZALES 

0  ICY 

ATTN 

D.  MCDANIEL 

Stewart  radiance  laboratory 

.TAH  STATE  UNIVERSITY 
*  jE  ANGELO  DRIVE 
sEDFGRD,  MA  01730 
01CY  ATTN  J.  ULWICK 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
LEDFORD,  HA  01730 

O’CY  ATTN  W.P.  BOQUIST 

TRW  DEFENSE  &  SPACE  SYS  GROUP 
CSS  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
0 ICY  ATTN’  R.  K.  PLEBUCH 
0 ICY  ATTN  S.  ALTSCHULER 
01CY  ATTN  D.  DEE 

Vi  SI  DYNE 

SOUTH  BEDFORD  STREET 
U.-E.INGION,  MASS  01803 
■-2CY  ATTN  K.  REIDY 
01 CY  ATTN  J.  CARPENTER 
C1CY  ATTN  C.  HUMPHREY 
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IONOSPHERIC  MODELING  DISTRIBUTION  LIST 
(UNCLASSIFIED  ONLY) 

PLEASE  DISTRIBUTE  ONE  COPY  TO  EACH  OF  THE  FOLLOWING  PEOPLE: 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 
DR.  P.  MANGE  -  CODE  A 101 
DR.  R.  MEIER  -  CODE  4141 
DR.  E.  SZUSZCZEWICZ  -  CODE  4187 
DR.  J.  GOODMAN  -  CODE  4180 
Dr.  R.  RODRIGUEZ  -  CODE  4187 

A.F.  GEOPHYSICS  LABORATORY 
L.G.  HANSCOM  FIELD 
BEDFORD,  MA  01730 
DR.  T.  ELKINS 
DR.  W.  SWIDER 
MRS.  R.  SAGALYN 
DR.  J.M.  FORBES 
DR.  T.J.  KENESHEA 
DR.  J.  AARONS 
DR.  H.  CARLSON 
DR.  J.  JASPF.RSE 

CORNELL  UNIVERSITY 
ITHACA,  NY  14850 
DR.  W.E.  SWARTZ 
DR.  R.  SUDAN 
DR.  D.  FARLFY 
DR.  M.  KELLEY 

HARVARD  UNIVERSITY 
HARVARD  SOUARF. 

CAMBRIDGE,  MA  02138 
DR.  M.B.  McELROY 
DR.  R.  LINDZEN 

INSTITUTE  FOR  DEFENSE  ANALYSIS 
400  ARMY/ NAVY  DRIVE 
ARLINGTON,  VA  22202 
DR.  E.  BAUER 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
PLASMA  FUSION  CENTER 
LIBRARY,  NW1 6-262 
CAMBRIDGE,  MA  02139 

NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREFNBELT,  MD  20771 
DR.  S.  CHANDRA 
DR.  K.  MAEDA 
DR.  R.F.  BENSON 


NATIONAL  TECHNICAL  INFORMATION  CENTER 
CAMERON  STATION 
ALEXANDRIA,  VA  22314 
I 2CY  ATTN  TC 

COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 
DR.  T.  CZUBA 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO.  CA  92152 

MR.  R.  ROSE  -  CODE  5321 

NOAA 

DIRECTOR  OF  SPACE  AND  ENVIRONMENTAL 
LABORATORY 
ROULnER,  CO  80302 
DR.  A.  GLENN  JEAN 
DR.  G.W.  ADAMS 
DR.  D.N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.  F.  DONNELLY 

OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  OUINCY  STREET 
ARLINGTON,  VA  22217 
DR.  H.  MULLANEY 

PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PA  16802 
DR.  J.S.  NISBET 
DR.  P.R.  ROHRBAUGH 
DR.  L.A.  CARPENTER 
DR.  M.  LEE 
DR.  R.  DIVANY 
DR.  P.  BENNETT 
DR.  F.  KLEVANS 

PRINCETON  UNIVERSITY 
PLASMA  PHYSICS  LABORATORY 
PRINCETON,  NJ  08540 
DR.  F.  PERKINS 
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SCIENCE  APPLICATIONS,  INC, 

1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 
OR.  O.A.  HAMLIN 
DR.  L.  LINSON 
DR.  E.  PRIEMAN 

STANFORD  UNIVERSITY 
STANFORD,  CA  94305 
DR.  P.M.  BANKS 

U.S.  ARMY  ABERDEEN  RESEARCH 
AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  MD 

DR.  J.  HEIMERL 

UNIVERSITY  OF  CALIFORNIA, 
BERKELEY 

BERKELEY,  CA  94720 
DR.  M.  HUDSON 

UNIVERSITY  OF  CALIFORNIA 
LOS  ALAMOS  SCIENTIFIC  LABORATORY 
J-10,  MS-664 
LOS  ALAMOS,  NM  B7545 
M.  PONGRATZ 
D.  SIMONS 
G.  BARAS CH 
l.  DUNCAN 
P.  BERNHARDT 

UNIVERSITY  OF  CALIFORNIA, 

LOS  ANGELES 
405  HILLGARD  AVENUE 
LOS  ANGELES,  CA  90024 
DR.  F.V.  CORONITI 
DR.  C.  KENNEL 
DR.  A.Y.  WONG 

UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20740 
DR.  K.  PAPADOPOULOS 
DR.  E.  OTT 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 
DR.  N.  ZABUSKY 
DR.  M.  BIONDI 
DR.  E.  OVERMAN 


UTAH  STATE  UNIVERSITY 
4TH  AND  RTH  STREETS 
LOGAN,  UTAH  S4322 
DR.  R.  HARRIS 
DR.  K.  BAKER 
DR.  R.  S CHUNK 


